Infinite Spatial Domains and the Fourier Transform 

18.303 Linear Partial Differential Equations 



Matthew J. Hancock 
Fall 2006 



[Dec 5, 2005] 

Consider the heat equation on an "infinite rod" 



Ut = nu^x, —oo < X < oo, t > 
u{x,0) — f {x) , — oo < X < oo. 

Since there is no boundary, we don't have so-caUed boundary conditions. However, 
we assume \u{x,t)\ and remain bounded as \x\ — > oo (we wiU make the re- 

quirements more precise soon). We proceed by separation of variables, 

u{x,t) ^X{x)T{t) 

On substitution into the PDE, we obtain, as before, 

r _X" _ 

and A is constant since the Lh.s. depends only on t and the middle only on x. The 
individual problems for X and T are 

X" + XX = (no BCs) 
T' = -kXT 



Solving for T gives 



T (t) = ce 



—nXt 



Wc have not said yet whether A is negative, positive, or zero. Since we the energy 
in the rod will stay the same, and will dissipate, we expect \u {x, t)| < oo for all times 
t > 0. If A < 0, T (t) — oo as t — s> oo and X {x) blows up too, and hence A must be 
non-negative. If A = 0, then T {t) = X = Ax + B, and A = to avoid a blowup. 
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In fact, since we assumed — > as — > oo, we must have ^4 = 0, so that X — 0. 
Thus A = leads to a trivial solution. It turns out that we can relax the conditions 
on / (x) to allow A = 0. 
For A > 0, we have 

X = Acos (yXx^ + B sin ^a/Ax^ 

It is customary to let A = a;^ (where a; > 0), so that for each A = cu^ > 0, we have a 
solution 

u {x, t; iu) = {A (iu) cos (lux) + B (u) sin (ujx)) e"^^*^* 

Since G M, we must integrate, instead of summing, over all possible uj, to obtain 
the full solution: 

poo 

u{x,t)^ / {A (cu) COS (lux) + B{uj) sin {ujx))e'''^^'^''duj (1) 
Jo 

Imposing the initial condition u {x, 0) = / (x) gives 

poo 

f (x) = / {A (cu) COS (cux) + B (cu) sin (oux)) duj (2) 
Jo 

The concept of orthogonality for a discrete set of eigen-functions can be generalized 
provided / (x) satisfies the following properties, 

/ (x) is piecewise smooth on every interval [a, h] of the real line 

1/ {x) \ dx < oo. 

The second condition ensures / {x) decreases fast enough for large If these prop- 
erties of / {x) hold, then A {uj) and B {uj) are given by 

A{uj) — — 1 f (s) cos (us) ds, B{uj) — — I f (s) sin (us) ds. 

TT J-oo TT J_^ 

Thus, 

-j^ /*cxj / noo 

u{x,t) = -J e""^'"* (cos (ux) J f (s) cos (us) ds (3) 

+ sin (ujx) / / (s) sin (cos) ds I duj 

J — oo / 

oo poo 

/ e~'^ (•s) (cos (ujx) cos (us) + sin (ux) sin (ujs)) dsduj 

J —oo 

/oo 
e"'''"*/ {s) cos (a; {x - s)) dsdw 
-oo 

e-'"''^* cos {u {x-s)) du?j f (s) ds 
K{s,x,t)f{s)ds (4) 



oo 

oo POO 



oo 

oo 
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where 



1 f°° 2 

K {s,x,t) = - e cos {cu {x - s)) dcu 



(5) 



is called the Heat Kernel. 

The integral for K can be calculated by defining z = ujy/Kt, so that 



K {s, X, t) — 



Let 



TTV K,t 



I{h) 



2 

cos {hz) dz, b 



e ^ cos (bz) dz 



X — s 



(6) 



You can find / {b) from Tables, e.g. Table of Integrals, Series & Products by Grad- 
shteyn & Ryzhik. But here's how to calculate it: first, differentiate in b, 



I' (b) = (~^~''^) ^^^) 



Integrating by parts with U — sin 62;, dV — —ze ^ , gives 

e 



I'{b) 



sin (bz) 







- / e cos (bz) dz 



Thus 

Integrating yields 
and from (6), 

Therefore, 



7 (6) = 7 (0) e-^'/^ 



7(0) 



e ^ = 



Substituting 7 (6) into (5) gives the Heat Kernel 



K (s, X, t) = 
and temperature, from (4), 



[X — s] 



exp 



Ant 



/oo i*c 
K (s, X, t) f (s) ds = 
-00 ^— c 



: exp 



— sj 

4:Kt 



ds. 



(7) 



(8) 



Since K {s,x,t) decays exponentially as , |s| — ^ 00 , then a sufficient condition for 
a solution is that |/ is bounded above by some constant as \s\ — > 00. 
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1 Hot spot on an infinite rod 



As an example, consider the initial hot spot 



uq/c, \x\ < e/2 
0, \x\ > e/2 



where £ > 0. The temperature distribution is 



u{x,t) — 



exp 



Mo 



e/2 



e/2 



exp 



[x — sY 

4:Kt 

(x — sY 
Ant 



ds 



ds 



To evaluate this, we note that the error function is defined as 



erf (x 

Let z ^ {x - s) /VAkI. Then 

f(a;-£/2)/V4Kt 



)=-r 



e dx 



u (x, t) 



Up r 



Uq 



£\/lT 
Uq 



exp (— dz 

(x+£/2)/V4/tt 

j-{x-e/2)/^/iKt 

exp (— dz + exp {—z^) dz 

{x+e/2)/V^ Jo 



e\/7r 



(x+£/2)/V4ret 



exp {—z^) dz 



ix-e/2)/V4M 



exp {—z^) dz 



2e 



erf I erf ^^'-^/^ 



Plots are given below in Figures 1 to 4. The solution u (x, t) is just the sum of 
shifted erfx functions, and hence is straightforward to plot once you know what erf(,T) 
looks like (see Figures 1, 2). The error function is 2/y/n times the area under e^*^ 
from s = to s = X (Figure 1). Thus as s — > oo, the area is -^/7r/2 since erf (oo) = 1. 
To obtain Figure 3, note that 



and 



Uq 



exp 



jx + €/2Y 

4:Kt 



exp 



4:Kt 



Ut < 0, \x\ < e/2 



4 
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Figure 1: Plots of e and eri{x). 
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Figure 2: Plot of t) for t = and Kt = 1, with s = 1. 

For I a; I > s/2, u increases and then decreases. Thus Ut = Q when 

x + e/2 ex 



In 



X - e/2 2Kt 



and hence when 



ex 



2k In 



x+e/2 • 
x-e/2 



2 Fourier Transform 

Ref: Guenther & Lee, §5.4 and for background, see §3.4 (in particular p. 78 and 79) 
Rcf: Myint-U & Debnath §11.1-11.4 
[Dec 7, 2005] 
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Figure 3: Plot of u{x, t) for x — (solid), x = ±1/4 (dash) and x = ±3/4 (dot-dash), 
for e — 1. 
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Figure 4: Plot of the level curves of u{x,t), i.e. the curves on which u{x,t) = const. 
Numbers indicate the value of u{x,t)/uo on the level curve. Here e = 1. 
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We now introduce the Fourier Transform and show how it is related to the solution 
of the Heat Problem on an infinite domain. Consider a function g {x) defined on 
— oo < X < oo that satisfies the properties listed above, namely, 

g {x) is piecewise smooth on every interval [a, h] of the real line (9) 

/oo 
\g{x)\dx < oo. (10) 
■oo 

The Fourier Transform G [uj) of the function g [x) is defined as 

G{u;) = ^J^g{x)e'-^dx (11) 
The Inverse Fourier Transform g (x) of a function G (a;) is given by 

/oo 
G (co) e-^^dw (12) 
-oo 

2.1 Fourier transform and the solution to the heat equation 

Ref: Myint-U & Debnath example 11.3.1 

To relate the solution of the Heat Problem on an infinite domain — oo < a; < oo 
to the Fourier Transform, we must make some manipulations to our solution. In 
particular, Eq. (1) becomes 



u{x,t) = I {A{uj)cos{ujx) + B{uj)sva.{ujx))e '^^djjj 
Jo 

POO POO 

= / Fi (cu) e*^^e-'^'''*(icu + / F2 (cu) e-'^'^e-'^'^'du, (13) 
Jo Jo 

where 

F, (a;) = 1{A (a;) - iB (a;)) , F, (a;) = ^ (^ (a;) + iB (a;)) . 

We now allow 00 to take all real values, i.e., —00 < uj < 00, and continue to manipulate 
the solution u {x,t). Making the change of variable a; — > —a; in the first integral of 
(13) gives 



/O /-oo 
Fi{-uj)e-^''e-'''^'du+ / F2 (a;) ^ 
00 Jo 

/oo 
F (u;) e-^'^e-'^'^'du;, (14) 
-00 



where 

F2(w), uj>0, 



F{u;) 



Fi (-tu) , a; < 0. 



Furthermore, ai t — 0, we have 



oo 



u{x,0)^f{x) 



(15) 



— oo 



Comparing Eq. (15) to (12) shows that / (x) is the Inverse Fourier Transform of F (uj), 
or equivalently, F [uj) is the Fourier Transform of the initial temperature distribution 
/ {x). From (11), we have 



Therefore, given an initial temperature distribution / (x) that satisfies conditions 
(9) and (10), we find the Fourier Transform F (u) of / (a;), and the temperature 
distribution at each point in time is then given by (14). The integral in (16) can 
generally be found in Fourier Transform tables. 

2.2 Fourier transform of the heat equation 

Ref: Myint-U & Debnath §11.2 

We now define the Fourier Transform (FT) of a function u {x, t) as an operator: 



Thus, the FT T maps a function of {x,t) to a function of {u!,t). To transform the 
Heat Equation, we must consider how the FT maps derivatives. Note that 




(16) 




(17) 




Also, integration by parts and the fact that u {x, t) 
calculate the FT of du/dx, 



as b 



oo allow us to 




iujU (a;, t) . 
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Thus, the FT of an x-derivative of a function is mapped to —iu times the FT of the 
function. Hence 



du 
dx 



-iu)^ !F[u\ — —u^U {u!,t) 



Thus the FT of the Heat Equation ut = kUxx is 

d 



dt 



Hence, the Fourier Transform maps the heat equation, a PDE, to a first order ODE! 
Integrating in time gives 

U{u;,t)^C{u;)e-^'' (18) 

where C (cu) is an arbitrary function, due to partial integration with respect to time. 
Setting t = in (17) and (18) gives 

-j poo 1 poo 

C (a;) = C7 (a;, 0) = — / u {x, 0) e^'^dx = — / f (x) e^'^dx. 
27r i_oo 27r 

Wc could substitute C (u) back into (18) and use the Inverse Fourier Transform to 
obtain u{x,t). The solution would be the same as that for separation of variables 
(recall that uniqueness holds). However, u {x,t) can be obtained almost immediately 
using a result called the Convolution Theorem. 

One important FT is that of the Gaussian distribution. 



Thus, the FT of a Gaussian is a Gaussian. In particular, the IFT of a Gaussian is 



2.3 Convolution Theorem 



(19) 



Ref: Myint-U & Debnath §11.3 

Convolution Theorem Suppose that F (uj) and G (a;) are the Fourier Trans- 
forms of / (x) and g (x) , then 

poo 

2^ 



/oo 
F{u;) 
■oo 



/ {x) e^'^^rfx, G{uj) = — g {x) e^'^'^dx, 

2n 



/oo 
G (u) e-^'^duj. 
■oo 
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Let H{uj)^F {uj) G {uj). The Inverse Fourier Transform (IFT) of H {uj) is 

^(^) = ^y_ f{s)g{x-s)ds=—J^ g{s)f{x-s)ds (20) 

The first integral in (20) is called the convolution of / (x) and g (x). In other words, 
the IFT of the product of two FTs is the l/27r times the convolution of the two 
functions. 

Proof: [optional] The IFT of H (ou) is 

roo poo 

h{x)= H (a;) e-'^^'dcu = F{cu)G {u) e-'^^^dw 



Substituting the IFT of F (cu) gives 



h{x) = I [ f{s) e'^^'ds ) G {^) e-'""da; 



27r 



1 

2^ 



oo /"OO 



/ (s) G {u) e-^''e''''dsdcu 



oo ^ — oo 

oo 
oo 



Assuming we can interchange the order of integration (we can provided \f {s) \ ds 
and \G (a;)| dcu are finite and / (s), G {lu) are smooth) we have 

1 roo / /"OO \ 

hi^) = ^J fis)U G (uj) e-'^^^-'^duj ds 

Note that 

/oo 
-oo 



and hence 



1 /"^ 

^ (a:^) = ^ y f{s)g{x- s) ds. 



This is the first integral in (20). To obtain the second, we make the transformation 
w = a; — s, so the integral becomes 



1 f°° 

h{x) = — g{w)f{x- w) dw. 

27r J_oo 



Since w is a dummy variable, we can replace it with s to obtain the second integral 
in (20). ■ 

To apply the convolution theorem to the Heat Equation, we note that the FT of 
the solution u{x,t) is 

U {u, t) = C {u) e-^^* 

Since the IFT of C (lu) is / (x) and the IFT of e""^'* is, by (19), ^tt/ (fi;t)e-^'/(^'^*), 
then, by the Convolution Theorem, 

which agrees with the solution we found using separation of variables, Eq. (8). 
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2.4 Heat equation on a semi- infinite domain 

Ref: Guenther & Lcc p. 171, Myint-U & Debnath example 11.2.1 
Consider the heat equation on a semi- infinite domain < a; < oo, 

Ut — KUxx, < X < oo 
u{x,0) — f (x) , < X < oo 

and f (x), u{x, t) approach zero fast enough as \x\ — > oo so that the integral 

\u {x, t) \ dx 



is finite. Also, a boundary condition is imposed a,t x — 0. Either a Type I BC (fixed 
temperature) , 

u{0,t)^0, t>0, (21) 

or Type II BC (insulated), 

Ux{0,t) = 0, t>0. (22) 

To solve this problem, we recall the temperature distribution on the infinite do- 
main — oo < X < oo due to an initial temperature / (x) was given by 



/oo ^ 
K (s, X, t) f (s) ds 
-oo 



where the Heat Kernel is defined as 

K (s, X, t) — , exp 

Note that at x = 0, the Heat Kernel is even in s, i.e. K (s, Q,t) — K (— s, 0, t). Also, 



/oo ^ 
K{s,0,t)f{s) ds 
-oo 



Thus u (0, t) = if / (s) is odd. Therefore, the solution to the Heat Problem on the 
semi-infinite domain < a; < oo with zero temperature ai x = {u (0, t) = 0) is 



/oo ^ 
K (s, X, t) f (s) ds 
-oo 



where / (s) is the odd extension of / (s), i.e. 



f{s) s>0 

f{s) = { s = 

-f{-s) s<0 
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Similarly, note that the x-derivative of the Heat Kernel is odd at s = 0, 



and hence the solution to the Heat Problem on the semi-infinite domain < a; < oo 
with an insulated BC at a; = {ux (0, t) = 0) is 

/"OO ^ 

u{x,t)= / K (s, x,t) f (s) ds 



where / (s) is the even extension of / (s), i.e. 

f{s) s>0 
f{s) = { s = 
f{-s) s<0 

3 Fourier Transform solution to Laplace's Equa- 
tion 

Ref: Guenthcr & Lee §8.2, problem 5, Myint-U & Dcbnath example 11.2.1 

Suppose the temperature of an infinite wall is kept at / (x), for — oo < x < oo. 
Find the steady-state temperature in the region adjoining the wall, y > 0. The 
steady-state temperature satisfies Laplace's equation, 

V'^ue — 0, — oo < X < oo, y > 0. 

The BCs are 

ue — f {x) , — oo < x < oo, y — 0, 
lim Ue {x, y) — 0, lim ue {x, y) — 0. 

y^oo |a;|^oo 

We employ the Fourier transform in x, 

J'[g(x,y)]^^J^ g(x,y)e'-^dx 
We define UE{oj,y) — T \u,e {x-, v)\- As before, we have 

[uex^ = [mb] ^ -lo^Ue (o;, y) , T [uEy-^ = {ue\ = ("^^ V) ■ 

Hence Laplace's equation for the steady-state temp ue {x,y) becomes 

^C/b {uj, y) - uP'Ue {oj, y) = 
12 



Solving the ODE and being careful about the fact that u can be positive or negative, 
we have 

Ue (o;, y) = ci {u) e-l'^l^ + C2 {u) e^^^^ 

where c\ (cu), C2 (tu) are arbitrary functions. Since the temperature must vanish as 
y — > oo, we must have C2 (cu) = 0. Thus 



Imposing the BC at y = gives 



(23) 



ci (a;) = Ue (c^, 0) = {ue (x, 0)] 
Note that the IFT of e-l'^l^ is 



2 /""^ 
2^ 



Me (x, 0) e^'^dx 



2tt 



f {x) e^'^dx. 



/oo pu poo 

-00 J —CO J 



+ 



-00 

1 



(7/ + ix) 
% 

y — ix y + ix x'^ + y'^ 



y — ix 
1 

— + 



Therefore, applying the Convolution Theorem to (23) with T ^ [c\ (a;)] = / {x) and 

jF-l ^g-|a;|j/j ^ 2y/ (x2 + y"^) givCS 



1 /""^ 



2^ 

{x - s)^ + 



-ds. 
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